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Abstract 



In two-dimensional lattice fermion model a determinant representation for the two- 
point correlation function of the twist field in the disorder phase is obtained. This field 
is defined by twisted boundary conditions for lattice fermion field. The large distance 
asymptotics of the correlation function is calculated at the critical point and in the scal- 
ing region. The result is compared with the vacuum expectation values of exponential 
fields in the sine-Gordon model conjectured by S.Lukyanov and A.Zamolodchikov. 



During the last time the progress has been made in calculation of the long distance 
asymptotics of a correlation functions of local fields in some integrable two-dimensional 
quantum field theories Jl], |2|. Usually in these theories two-point correlation functions of 
local fields can be represented as an infinite series of the form-factor contributions, which are 
calculated using a method of the form-factor bootstrap || or the angular quantization 0. [5| . 
In [[], |7| a summation method of the form-factor decomposition of the correlation function 
was developed. This method allows to obtain a closed expression for the correlation function 
through Fredholm determinant of the integral operator and to do analysis of asymptotic 
behaviour of one. 

In this paper we analyse a large distance behaviour of the two-point correlation function 
of the twist field /^(r) in the disorder phase of two-dimensional lattice fermion model. At 
v — ~ this field coincides with the disorder field /i(r) and the correlation function of these 
fields satisfies the following relation || |J : 



1. Introduction 
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where // z (r) is a disorder field in the Ising model. 

Using the functional integral method in two-dimensional lattice fermion model, we obtain 
a determinant representation of the correlation function. Calculating the asymptotics of the 
correlation function for r — > oo, we find the vacuum expectation values (fj, u (r)). Here index 
v (0 < v < 1) denotes twisted boundary conditions for lattice fermion field along the 
correlation line. In scaling limit (massive fermion quantum field theory) this field is massive 
analog of the twist field in two dimensional conformal field theory with central charge c = 1 

In the free fermion point of the sine-Gordon model the vacuum expectation value 
< yu(r) > is connected by simple relation with the vacuum expectation value of the ex- 
ponential field exp(iu(f)) for v = ~. The sine-Gordon model is described by the Euclidean 
action 

1 f 

S = — dxdr (dk4>dk4> + /icos(/30) 

4:71 J \ 

The free fermion point occurs for (3 = 1. 

In this point the exponential field can be defined by the braiding relations with the free 
fermion field 01: 



jv(f>(T,x) 



As it was shown in 
following relations 

From here one gets 



^(r, y) = e 27riu ^{r, y)e iv ^ T ' x \ if y > x. (1.2) 
F2| the correlation functions of the disorder fields satisfy the 



(//(0)//(r)) = (cos 10(0) cosi0(r)>. 



(//(0)//(r)) >= i(e^°) e*3*W) + ]^{S m e^M). 



;i.3) 



;i.4) 



Using this equation, we can obtain a relation between the vacuum expectation values of the 
disorder and exponential fields. 

Note that for calculation of the left and right hand sides of ( |1.4| ) it is nessesary to 
regularize explicit expressions for the correlation functions that it leads to appearence of 
the corresponding wave renormalization constants. To be rid of them at calculation of the 
expectation values it is convenient to normalize the correlation functions on behaviour for 
r — > 0. In this limit we obtain from ( |1.4| ) 



(MOV(r)) = he 



;i.5) 



where we have used the following asymptotics for the correlation function of the exponential 
fields for r — > lT3l 



Jv(t>(0) „iv'cj>(r) 



> = <e 



i(v+v')ct>\ r 2w' 



1.6) 



Then for a ratio of (|1.4j ) in the limit r — > oo and (|1.5|) one gets 



As it will shown below one can choose such normalization of the disorder fields that asymp- 
totic behaviour of the correlation functions in denominators of the left and right hand sides 
of ( |1.7| ) coincides. In this case we obtain from ( |1.7| ) (n) = y / 2(exp(i|0)). In the paper |T3 | 
(also see ]14|) it was conjectured the following expression for the vacuum expectation values 
of the exponential fields (\u\ < 1) 



{e iu<f> ) 



exp 




This expression it is easy to obtain by means of results of the paper [plj] . Here the asymptotic 
behaviour of the correlation function of the exponential fields for r — > was calculated. 
Comparing our result ( |4.20| ) and (|1.8j), we obtain 



2 

(fi u ) = A v (exp(«V0)}. 



;i.9) 



Emphasize that this relation is correct for such normalization of the fields that the correlation 
functions of the twist and exponential fields have the same asymptotics for r —>■ 



</^(0)Mr)> 



1 
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In section 2 the determinant representation of the two-point correlation function of the 
twist field is derived. In section 3 asymptotic behaviour of the correlation function are 
considered and new determinant representation through the matrix of the Toeplitz type 
are obtained in the limit r — > oo. In section 4, using this representation, asymptotics of 
the correlation function is calculated. In Appendices some relations used in this paper are 
derived. 



2. Functional integral representation 
of the correlation function 



It is known Jl5|, |I6| , that the correlation function (/i / (0)/i / (r)) in the two dimensional 
Ising model can be represented in the form of the functional integral for the lattice Majorana 
fermion theory with antiperiodic boundary conditions for the fermion field with the exception 
of the line [0,r — 1], where the fermion field has periodic boundary conditions. 

The same way let us define the two-point correlation function of the field (i u (r) in the 
lattice Dirac fermion theory through a ratio of the following functional integrals 



with the action 

Srty] = (W) = E^(p)iVVV), 

p,p' 



(2.2) 



where ip{p) is a complex two component grassmann field; the coodinates of the lattice sites 
p = (x, y) run through the values x = 1, n x , y = 1, n y ; D is the lattice Dirac operator 



n_ u v \ _ f l-tV x l-tV y \ ,, 

Here V x , denote the shift operators along the X and Y axes: 

VxV(p) = i>(p + x), V y V»(p) =ip(p + v), 

where x, y are unit vectors. Index "r" denotes that in action S r {ip] operator V y has twisted 
boundary conditions along the line [0, r — 1] (r — 1 = (r — 1, 0)): 

V y ip(x, n y ) = -e 2 ™^p(x, 0), x = 0, 1, . . . , r - 1. 

In denominator ([2.1|) the action S[ip] coincides with S r [if>] by the form, but Vj, satisfies 
antiperiodic boundary conditions along all boundary. 

It is not hard to check that the operator D satisfies the relation D^ 1 = D T ■ K, where K 
is the lattice Klein-Gordon operator 

K = detD = uu T + vv T = 2(1 - t) 2 - t[(V x - 2 + V-*) + (V„ - 2 + V_„)]. 

In the "naive" continuum limit (the lattice constant a — > 0) K = a 2 t(m 2 — didi), where 
di = (Vj — l)/a and ma/y/2 = £~2 — tz. From here the critical point of the model is 
t c = 1, and the scaling region is determined by relation \t c — t\ = ma/y/2. In momentum 
representation the operator K is diagonal and the function 



K(p) = 2(1 - tf + At (sin 2 £p + sin 2 |) 



has unique minimum at p x = p y = in the Brillouin zone and unlike the usual lattice Dirac 
operator we have not problem with the fermion dubling. 

In order to integrate over the fermion fields in the numerator ( |2.1| ), in the action S r [i/j] it 
is convenient to go to operator V y with antiperiodic boundary conditions along all boundary. 
After that we get additional term 5S r [i()] in the action S r [i/)], which contains a lattice defect: 

Sr[ll>]=S[l/>)+6S r [ll>], 

where 

5S r ty] =tJ2{-&\x,n y )^ 2 (x,0)+C^ 2 (x,0)^(x,n y )] = (^P T EP^). 

When 

<M0K(r» = = ■ (2.4, 

Here the projective operator P selects lattice sites on the line [0, r — 1] 

Pl,x = $l,x5 ny ,y ( J y J, / = 0, 1, . . . ,r - 1; 



and 

-£7» 



£*j(r) g 



where £ = 1 — e 2mu and 7^ is a unit matrix of the size r x r. 

For calculation of the functional integrals in (|2.4j) it is convenient to represent e' 5Sr ' 
through the integral over a auxiliary field 



e sSrW = |s| y exp^S" 1 *) + (xA/>) + (V^ T x)]- (2-5) 
Using(2.5) and integrating over ip(p) in (|2.4j) , one gets 

(Mo)MO> = 1=1 /d|xx]exp[(xs- 1 x) + (x^~ 1 ^ r x)] = 



Id Ic— 1 i Dn-loT 



PD- A P J | = |G|, (2.6) 

where \G\ denotes the determinat of a block matrix of the dimension (2 x 2) with blocks of 
the dimension (r x r). The matrix G has the form 



, i ^ sin ttz/ -gf , sin vrz/ + zi£], cos ttz/ , 

U ~~ 12 • , tM / 11 \T • ' v 

Sin 7TZ/ + Zi^ ^, COS 7TZ/ U7i,a:') sm ^ v 



where 



ii 1 f d 2 p 

9x ' x ' (2vr) 2 J K{pf 



Jp*(x-x')2tu*(p x ), u{p x ) = 1 - te ip * 



— 7T 

Let us show that the correlation function ( |2.6j ) for v = ~ is equal to square of the correlation 
function of the disorder field in the Ising model. For this it is necessary to do a similarity 
transformation of the matrix G in fl2.6| ) by means of the unitarity matrix Q 

Q = eK 75= (° "*)' Ctg ^ = t ' (2 ' 8) 

As result we obtain 

G{r) = (/i^(0)^(r)) = \Q-G-Q- 1 \ = \ cos 2 vrz/ + sin 2 vrz/V • V T \ } (2.9) 

where 



2h j z \ (1 — az — pz) 

I J2 1 = 1 ' 



(2.10) 



a = t— = , p = t— = . 2.11 



Here the matrix V x>x > exactly coincides with classical expression for the Toeplitz matrix |L7| . 
the determinant of which determines the correlation function (p 1 (§) p, 1 (r)) in the paramag- 
netic phase of the Ising model. It is obvious that ( |1.1| ) follows from Q2.9| ) for v = |. In 
further we include coefficient (2iri)~ 1 in the integration measure of contour integrals. 

Note that the transformation Q is connected with particular version of the lattice Dirac 
operator (2.3). In the "naive" continuum limit we obtain for the action (2.2) 

Sty] = J d 2 p^(p)(me^ - 7 iW(p), ( 2 - 12 ) 

where 

7 *=(J -l)' 7 *=(z ~o)' 

At the critical point t c = 1 the angle 9 in ( |2.8|) is equal to it/ 4 and the transformation Q is 
the 75-rotation of the Grassmann field: ip' = exp(ry 5 7r/8 )■?/;, after that the Dirac action takes 
the usual form 

Sty] = J d 2 P 4>(p)(m - 7 iW(p)- (2-13) 



3. Toeplitz determinant representation 
of the correlation function 



correlation function (|2.9| ) can be reduced to a calculation of the Toeplitz determinant: 



where 



In this section we show that in the scaling limit (ma — > 0) at r —>■ 00 a evaluation of the 

;ion of the Toeplitz < 

y(y) . v ( * u)T \ = |V M | 2 , 
(l-az)(l -I3Z- 1 )' 



G(r) 



cos ttu + sin nvV ■ V 



(3.1) 



dz 



— z -(*-*')vM(z), V {u \z) 



z=l 



(3.2) 



is a Toeplitz matrix. For v — ~ the kernel V^ v \z) coincides with the kernel V(z) ( |2.10| ). 

In Appendix A it is shown that the correlation function (|2.9| ) can be represented in the 
following form 



G(r) 



|1 — sin txvA\ , 



(3.3) 



where 



A XjX > = j 



dz 1 dz 2 (z 1 nz 2 r' v ^ , 

(1 - z 1 z 2 ) 



\Z1Z2\<1 

Here the integration contour passes between the cuts depicted on Fig. 1. 



z 




Fig.l. 



Using (3.3), we get 



afc /sin7rzA 



2k 



ilnG(r) = Spln(l - sin 2 7rzM) = ~E y 



(3.4) 



k=l 



where 



Ctfc = 7T bp/i = 7T 



'li n^[(i 



nL[i-W 



2 2Z+1 — Z\- 



For r — > oo the terms in the right-hand side of 
are exponentially small, so that 



ziZi + i)V(zi)} 

which contain the product (z/^+i 



(3.5) 



a k = n 



2k 



n^[(i-^ + i)^,)] 



7 ng 
i nEi(i- 



2fc 

n 



(a - - (3zi) 
(1 - a2|)(z f - (3) 



(3.6) 



At ma = (a = 1) the integral in ( |3.6| ) has logarithmic divergence on the upper limit 
(in a neighborhood of the "upper right cone" of the 2/c-dimensional hypercube with the 
coordinates z\ = 1). Let us isolate and evaluate this divergence: 



2k dz, 2k ( c 

n^T— n(i 

i=i 1 - z i z i+i i=i u 



a - zi \ 2 



- OLZl 



e/ma 



2k 



n 



(3.7) 



For derivation of this estimate we did the following substitutions of the integration variables: 
z\ = a — z' h z[ = (majz", z" = zi + 1. Isolating the logarithmic singularity in the last integral, 



e/ma 2fc 



we obtain 



Differentiating over ma the left- and right-hand sides of (|3.8|) , the coefficients b k can be 
expressed through the multiple integrals 



dz 

a k — I TT = —kbk\n(ma) + const. (3.8) 

J iJ[ Zl + Zi+i 



1 2k-l p 2k-2 -,-1 

h = 2 J Y[dz r (1 + z x ){l + z 2k _ l ) J] + z l+1 



i=i i=i 



(3.9) 



Using ( p.8| ) and ( |3.4j ), we get the asymptotic evaluation of the correlation function for r — > oo, 
ma — > 

ilnG(oo) = ln(ma) ^ 6 fc ( ) + const. (3.10) 



7T 



fc=l 

The integrals in the first three coefficients in (|3.9|) can be reduced to the table ones 

h = l, 6 2 = 2C(2), 6 3 =4C 2 (2) + 6C(4), (3.11) 
where ((h) is the Rieman ^-function 



k\ J 1 — x 

o 

Decomposing in (|3.10| ) the functions (sin7rz//7r) 2A; over powers of v and using the values ( |3.11| ) 
for the coefficients b k , it is not difficult to get 

sin 2 irv 2 4 

, sin 2 7tz/ , sin 4 nu 9 ^ . K . , 

6i o— + &2 — = i/ 2 + 0(z/ 6 ), (3.12) 

, sin 2 7tz/ , sin 4 7rz/ , sin 6 ixv 9 ^ . s , 
6i 2— + b 2 — + h — = v 2 + Of 8 , 

7T 7T 4 71"° 

i.e. every following term in the series (|3.11| ) annihilates more high powers of u 2k . 



Although the integrals (|3.9| ) is hardly calculated for multiplicity higher five, one can 
assume that in more high orders relations of the type (|3.12|) also take place. Then 

x-^ , / sin TTU\ 9 

X>(— ) =" 2 - < 313 > 

It means that the generation function for the coefficients b k has the form 

J oo 

— arcsin 2 (7r v / ^) = V &fc-2 fc . (3.14) 
77 fe! 



Then for the determinant 03.11) we obtain 



G(oo) = const • (ma 

(rao)^O 



(3.15) 



Note that ( |3.14 ) allows to find values of the integrals ( |3.9|) of arbitrary multiplicity. 

Let us prove the asymptotics ( |3.15| ). The following observation is key for the proof (it 
follows from the calculations in ( |3.7| )): the asymptotic evaluation (ma —>■ 0) of the integrals 
in ( |3.6|) is insensitive to the exponent in the expression for kernel of the matrix V X)X r. It 
allows to introduce the new matrix V^}., with the kernel (|3.2| ). Let us define the correlation 
function 

G w (r) = ■ V^ T \, (3.16) 

where the matrix is determinated in (|3.2| ). Using the results of Appendix A, (|3.16j ) can 
be represented as 



where the matrix has the form 

(„) / dz 1 dz 2 (z l ) x (z 2 ) 



A (l/) \ 2 , 



A y > 



(3.17) 



(3.18) 



(l- Zl z 2 ) V^(z 1 )V^(z 2 y 
The trace of the k-th power of the matrix is expressed through the 2/c-multiple integral 



s P [(i ( 



u)\k] 



sin txv 



71 



2k " 2k 



dzi 



2k 

\ (1 — ZiZi+i) ^ 



(a - zi)(l - (3zj) 
(1 - azi)(zi - (3) 



v k 



Ull - ( Zl z l+I y]. (3.19) 
1=1 



Here, unlike (|3.4f) , the factor (sin7r^) 2fc appeares on account of the exponent v at contraction 
of the integration contour. Taking into account ( |3.19|) for r —>■ oo, we get the expression for 
G^Hr) which is similar to (B .41) 



ilnG M (r) = Spln(l-i M 



£ 

fc=i 



» 



sin7ri/ 



7T 



2/c 



where 



(") 



2k 







da 



y (1 - ^ [(1 - - /3)J 



2k 

n 



nti-te+i)' 



(3.20) 



(3.21) 



!=1 



For r — > oo the terms in (|3.21|) containing the power of (ziZi + i) r disappear. Although the 
rest integrals in ( |3.21| ) diverge on the upper limit the difference between ( |3.6| ) and (|3.21| ) 



Aa 



2k 



071 



i (i-^i+i) h=i 



(1 - Qt^)(^j - /3) 



/J [(1- a*,) («,-/?) 



(3.22) 



is finite for a 



2k j 2k 

Aa^ = V2(§ - v) f 71 V E( X ~ z + const = const 

x l 6 i=i I 1 ~ Z i z i+V i=i 



Hence a ratio of the determinants ( |2.9| ) and ( |3.16| ) is also finite and it is not equal zero 



7 = 



G(oo) 



cos 2 7tz/ + sin 2 TTuV ■ V T \ 



oo 



T\ 



(ma)-»0 



const. 



(3.23) 



Therefore the singular over (ma) factor in G(oo) coincides with one in the determinant 
(|3.16|) . Since the matrix V^ v ' is a Toeplitz matrix we can calculate (|2.9| ) in limit r — > oo, 
(ma) — ■> 0. 

Note that we have considered the ratio ( |3.23|) in limit r = oo, (ma^ 



0. In the paper [13 



the correlation function of the exponential field was normalized on the correlation function 
in the massless fermion field theory. Therefore for comparison of our results with |13| it is 
necessary to evaluate the ratio (|3.23|) in limit (ma) = 0, r — > oo: in other words we must 
verify that the value of 7 does not depend on a sequence of the limits. 

It can show that at (ma) = the difference of the coefficients (|3.5| ) and (|3.21|) differs 
from ( p. 22 ) on some term decreasing with increase of r 



(y) 
a k — a k 



<t>k(r) 



For 4>k(r) we have the following asymptotic evaluation 

ln fc r 



k\ 



As a result we get at ma = and r — > 00 

G(r) 
GM(r) 



7 exp 



EMr) 



k=l 



(3.24) 



(3.25) 



Using ( |3.24| ), it is not hard to sum the series in the exponent of the right-hand side of ( |3.25| ) 
and to obtain the following estimate 



G(r) 
G<y)(r 



7 + 0(l/r). 



(3.26) 



4. Asymptotics of the correlation function 

In this section we calculate the large distance behaviour of the correlation function G(r). 
We use the determinant representation (|3.1| ) for it in this case. Since in ( ^lf) the matrix v} u J, 
is a matrix of the Toeplitz type it can apply the McCoy and Wu techniques [0, which they 
used for the asymptotic evaluation of the correlation function in the Ising model (v — |). 
This techniques requires some generalization for < v < 1. In Appendix B we obtained the 
following asymptotic expression for the determinant of the matrix V"W (the formula |B . 24|) : 

\n\V M \ r =Jr + l) j^lnV^ - j dz\nP(z)^-\nQ( Z - 1 ). (4.1) 



Substituting here the explicit expressions for V^)(z), P{z) and Q(z): 



VM(z)=P(z)Q{z- 1 ), P{z) 



1 — az 



Q(z) 



l-(3z 
1 — az 



P(z) 



(4.2) 



one gets 



j ^ i n = vU ( ^]n P{z) - j> y In P(^ 1 )^) = i/ (P(0) - P(0)) = 0, 



- dz\nP(z)-^-\nQ(z' 1 ) = v 2 j> dz\nP(z) 



a z — (3 



(lnP(a) -lnP(/3)) 



Thus 



P(a) 



(1 - a/3) 



(4.3) 



Note that this expresion is derived for finite values of (ma) and r. The derivation of 
( |4.1|) is also correct in the scaling region: (ma) — > 0, r — > oo, r(ma) = const. However 
for (ma) = 0, r = const the integrals in the right-hand side of (|4.1|) are divergent and the 
calculation of | V^^) | at (ma) = and r — > oo requires the special consideration. 

Note that the following ratio of determinants 



(1 - a/3) 2 



(ma)=0 



(4.4) 



is finite both in the scaling regime and at the critical point. Here we explicitly indicated a 
dependent of the matrix on the parameter /3 (the matrix V'^(O) is determined by the 
kernel Q at (3 = 0). 

The relation fl4.4| ) it is not hard to obtain by means of ( [4.1| ). Really, for the ratio of the 
determinants in the left-hand side of (|4.4|) we get 



In 



j> dz\nP(z)—\nQ(z- 1 ) 



d 



+ Idz [lnP(z) + z/ln(l - (5z)\ ^- lnQfz -1 ) - u\n(l - fiz' 1 ) 
J oz L 



-v In P{(3)+v In P(0) - v l m(l - + 

+z/ IdzMl -I3z)— lng^" 1 ). 
J oz 



(4.5) 



The last term in the right-hand side of (|4.5|) one integrates by parts 

vj>dz\n{\ ~ (3z)-^\nQ(z- 1 ) = v(3 f -^—]nQ(z) = u\nQ((3). 
As a result we have 



In 



.|VM(0)|. 



-i/ In P(/3) + z/ In Q(/3) - v 2 ln(l - /? 2 ) = In 



1-/3 2 
;i - a/3) 2 



For (ma) = (a = 1) matrix elements of the matrix 



dz 



13=0 



are calculated in the explicit form 



(„),„, sin 7^ 



(4.6) 



7T Z/ — X + X' 

Call to mind that x, x' take the values 0, . . . , r — 1. To simplify notations let introduce 

Ur = V£>(0), (4.7) 



where the lower index for the matrix (^4.6| ) indicates its dimension. Multiplying the matrix 
( |4.6| ) from the left by the triangular matrix L r and the right by the triangular matrix R r , 
we get the following relation 



-U r -! J ' 



(4.8) 



where 



/ 1 

1 V 



i i 



i 



o \ 





V 



i 

r-l 



f + 1 

V o o 



-— \ 

r-l \ 




(4.9) 



j — 1 ' 



... ^T-l 7 

r— 1 / 

The determinant of a product of these matrises it is not hard to calculate 

~ * (-l)*"- 1 shitty r(r + i/)r(r - i/) 

1 r r| ~ * YHf) • 

This allows to get from ( |4.8| ) and ( |4.10 ) the simple recursion relation for the determinants 

. T(r + v)Y(r — v) 



(4.10) 



which has the solution 



\Ui — i — \U r 



\u r \ = n 



P(r) 



Using the following representation for T-function 



(4.11) 



lnr(z) 



dt ( e z — e 



t V 1 - e-* 
o x 

the solution ( 4.111) can be represented in the form 

\n\UJ = - 



+ (z-l)e 



-t 



dt sh 2 — 

t l j sh 2 4 



(4.12) 



For r — > oo the integral in the right-hand side of ( [4.12| ) is divergent for the zero. Rewriting 
(|4.12 ) in the form 



\n\Ur 



—v lnr — / — 



dt fsh 2 4 







t V sh 



2 / 



2 - v\- 1 



v 



2 / o - \— 2 / 

we can isolate this divergence. Now the integral depending on r is well defined and we get 



lnlfL 



—v lnr — 



dt ( sh 



2vt 



t \ sh 



2 t 



2 -zAr* 



+ 0(l/r 2 ). 



(4.13) 



Taking into account flOp, (fO|), ( |4.4|) , ( |4.7|) and ( |4.13| ), G(r) one finds the following asymp- 
totic evaluations for the correlation function 



G(r) 



r— >oo, a<l 



[l- a 2 ){l-(3 2 
(1 - a/5) 2 



2i^ 



(2^ 



ma 



,2^ 



G(r) 



r— >oo, a=l 





2u 2 


CO , , 


exp 


X 


_r(l-a/?) 2 _ 



l 2 ^ 
1 2 



2 - z/ 2 e-* 



2 U e 



2 -2pM 



.2y 2 



(4.14) 
(4.15) 



where 



t \ sh 



2 t 



2(ma) = 1 — a 2 

a— ►! 



Using these asymptotics one can assume that in the scaling limit (r — > oo, ma — > 0, mar 
s = const) the correlation function G(r) has the scaling form 



F(s, v) = \imr 2lj2 G(r) = 2 3 » 2 S 2u2 f(s, u), 
where the function f(s, v) has the asymptotics 



1 for s — >• oo, 

2 ~2» 2 e -2p(u) s -2u 2 for s 



0. 



(4.16) 



(4.17) 



This scaling behaviour follows from the arquments: 

1) the two-point scaling function F(s,u) have to interpolate between the behaviour at 
large distance away from the critical temperature (r — > oo, ma — > 0, so that s — > oo) and 
the behaviour at large distance at the critical point (r — > oo, ma — > 0, so that s — > 0), that 
is the asymptotics ( |4.17[ ) have to result in asymptotics ( |4.14| ), ( [4.15| ). 

2) for v — | the scaling function ( 4.16|) and the asymptotics (|4.17 ) have to coincide with 
square of the scaling function and the asymptotics for the two-point correlation function of 



the disorder field in the Ising model [18], |19| PQ] , 

It is convenient to define the ratio h(Rm), which does not depend on a normalization 
and the lattice cutoff a, 



h(Rm) 



G(R)\ 



R— >oo 



F(s,v) 



G(R)\ 



F(s,v) 



U->o 



(2mRy u exp 



2 

L 



dt ( sh 



2vt 



t V sh 



2 t 



z/ 2 e-* 



(4.18) 



where R = (ra) denotes a dimensional distance. 

To compare this asymptotics with the result of the paper [13] in the free fermion point 
we use the following normalization for the correlation function for R — > [TBI 



G(R) =R~ 2v \ 
Using this normalization, from the ratio ( f4.18| ) we obtain 



G(R) = U v f = (2m 



2u 2 



cxp 



2/ 7 



dt ( sh 



2ut 



sh 



2 f 



2 - v 2 e~ l 



(4.19) 



(4.20) 



5. Conclusion 

In this work a determinant representation for the two-point correlation function of the 
twist field is obtained in the lattice fermion field model. The large distance behaviour of 
this correlation function is calculated and the vacuum expectation values of the twist fields 
is found. This value differs from the vacuum expectation value of the exponential field by 
the constant. 

For explanation of this fact one can bring the following arguments. It is known that 
the quantum fields defined by the commutation relations (|1.2j) realize the operator solution 



of the isomodronomic deformation problem for the Dirac equation pl[. The correlation 
function of the twist fields fl2.1| ) one can interpret as the functional integral in the Dirac 
fermion theory with given monodromy properties for fermion fields. Then this correlation 
function can be connected with other solution of the isomodronomic deformation problem 
for the Dirac equation. In order to check this asssumption it is necessary to calculate the 
determinant (|2.9|) and to obtain a differential equation for the correlation function of the 
twist fields. 



Note that with our point of view in Appendix B of the paper [|TJ| instead of the vacuum 



expectation values of the exponential fields the vacuum expectation values of the twist fields 
was calculated but since the normalization on behaviour of the correlation function of expo- 
nential field for r — ► was used the right result (US) for the vacuum expectation values of 
the exponential fields was obtained. 

The authors thank N. Slavnov, S. Pakuliak, S. Lukyanov and M. Lashkevich for some 
useful discussions. One of us (V.Sh.) thanks A. Morozov for his hospitality at the Institute 
of Theoretical and Experimental Physics and Prof. G. von Gehlen for his hospitality at 
Physikalisches Institute d. Universitat Bonn. 

This work was perfomed with financial support from Ukrainian Found Fundamental 
Reseachers (project No. 2.5.1/051) and the INTAS program (Grant No. INTAS-97-1312). 

6. Appendix A 

In this Appendix we obtain the following representation for the determinant in ( [3.1|) for 

r — > oo 

g{r) = | cos 2 tcis + sin 2 iruC • C T \ = |1 - sin 2 nvA\ 2 (A.l) 



where for convenience we introduced the matrix 

"(1 -az)(l -fiz- 1 ) 



{l-az- 1 )^- (5z) 



(A.2) 



It has dimension (r + 1) x (r + 1) (x, y — 0, 1, . . . , r) and < fi < 1. For \i=\ this matrix 
coincides with the matrix V x , y in the determinant representation of the correlation function 
(|2.9|) . In QA.lp the matrix A has the form 

|«1«2|<1 



a; 



Let us represent the product C ■ C in the following form 

r \t\<i<\,\ z£(l - z/0 C(0 J z Z f\z\<\i\<x z£(l - z/0 C(0 " 

8(x - y) - A xy - B xy , (A. 3) 

where 

= _r d(zQz-^y C(z) r djzQz*? 1 
^ T\Z\<K\z\zZ{l-z/Z)C(Z)z-+l/J l-zt c(z)c(0 

In ( |A.3| ) we summed a geometric progression over x' and used C(z~ 1 ) = C~ 1 (z). Substituting 
( pOD in ( PP| ), we get 

</(r) = |1 -sin 2 vrz/(i + 5)| = |(l-s 2 i)-(l-s 2 5)-s 4 i-5|, (A.4) 

where s = sin7riA 

Note that matrix elements of the matrices A and S have the following asymptotic be- 
haviour 

A „ x+V D „ 2r—x—y 

A xv ~ a , jD X y ~ a . 
and therefore the product A ■ B is a exponential small matrix. Really, 

d( Zl ... ZA ) c(z 3 )c(z A ) (z 3 y+i - ( Z2 y +1 

(1 - Zl z 2 ){\ - z 3 z A ) C( Zl )C(z 2 ) {Zl> KH) z 3 (l - z 2 /z 3 ) ™ 

arlnrff ^ -C(z^)) ( i P^C(z)) (A.5) 



1 — az I \ J 1 — olz 



Taking into account ( |A.5 ) in ( |A.4| ), we obtain at r — > oo 

9 (r) = \(l-s 2 A)\.\(l-s 2 B)\. 



(A.6) 



Let us show that |(1 — s 2 A)\ — |(1 — s 2 B)\ for r — > oo. For this we write the matrix B in 
the form 

B = J-D-J, 

where 

/ ... 1 \ 
... 1 



J xy = 5(x + y-r) 



(J ■ J) xy = S(x-y), 



D 



V 1 ... J 

d(z£)z x £ y 



j-y 



l-z£ 



C(z)C(0- 



Using this representation, it is not hard to show that 

1(1-^)1 = 1(1-^)1 
The determinants in (|A.8|) it is convenient to decompose in the series 



(A.7) 



(A.8) 



lnjl - s 2 A\ = Sp In (l - s 2 A) = - £ — a 2fc , 



oo ^2k 



k=l 



oo ^2k 



In |1 - s 2 t)\ = Spin (l - s 2 b) = - -r- d 2k , 



where 



k=l 



a 2 k = Sp (A) , d 2 k = Sp (/)) 
For r — > oo the matrix A and l) can be represented in the factorized form 

A^F® ■ F {a) , D^F® ■ F {d) , F<$ = jdzz x+y C{z), F^ = j> dzz x+y C~\z), 



so that 



'dztC-^ztY 



i=i 

2k 



ZlZi + i 



d 2 k=(SpF^) 2k = fj[[-l 

Consider, for example, the case k — 1 

dz\dz 2 



( dziC(zi) 



i=i \ 1 ~ 



^2fc+l = Zl, 



(A.9) 



zi«a|<l (1 - ZiZ 2 ) 2 



C( Zl )C(z 2 ) 



dz\dz 2 C{z\) 



Z2^1/Z2 J\zx\<\Z3\ (Zl - Z 2 ) 2 C(Z 2 ) 



dz\dz 2 C(zi) f , _ . . d _ ] ( 

7 ^T^H + f dzC ^7T C ( 

|«a|<|*i| {zi - z 2 ) 2 C(z 2 ) J az 



(A.10) 



Such as in the right-hand side of ( |A.10| ) the second term is a integral over a closed contour 
from a derivative it is equal to zero 



j> dzC(z)^-C-\z) = - j dz^- In C(z) = 0. 



d 



In the first term in the right-hand side of ( |A. 10 ) we did the replacement z\ 



do 



dz\dz 2 



hz xZ2 \<\ (1 - z 1 z 2 ) 2 C(z 1 )C(z 2 ) 
For k > 1 it is not hard to get the recurent relation 

d 2 k = Q-lk + <^2fc-2 — «2fc-2- 



a 2 . 



(A.11J 



(A.12) 



For this in the integral representation ( |A.9| ) it is necessary to do the replacement zi — > 
After that we obtain new integration contours with > 1. Contracting sequen- 

tially these contours over the variables z\ so that |^^+i| < 1 and calculating residues in the 
poles zi = (zj+i) -1 , we get ( |A.12| ). 

Taking into account the "initial condition" ( |A.11|) , it is easy to check that a solution of 
the recurent relation ( |A.12| ) is a 2 k = d 2 k- From here we obtain ( |A.8| ). Then ( |A.8| ) and ( |A.6| ) 
lead to flOP . 

Note that putting \i = | in ( |A.2| ), we get ( |3.3| ) from ( |A.1| ), and putting v = h in ( |A.1| ) 
and /z = z/ in ( |A.2| ), we get (|3 . 1 7| ) from ( |A.1|) . 



7. Appendix B 

In this Appendix we obtain the asymptotic evaluation of \V^\ for r — > oo 

'(1-« 2 )(1-/? 2 )' 



It is convenient to introduce the following notations 



K« = B my = f -jz- x+ yB(z), B(z) = P(z)Q(z- 1 ), 
where x, y — 0, 1, r, the functions P(z) and is anylitical in the circle \z\ < 1 



P(z) 



az 



f3z 



Q(z) 



13 z 



az 



P(z) 



a > (3 < 0, 



(B.2) 



(B.3) 



a and (3 are determined in ( [2.11| ). 
Let us define 

dz 



B xy (fi) = j> -l z -x+y B »(z), B»{z) = P^Q^z- 1 ), 



(B.4) 



and 



f(fi) = hi\B(ji)\, /(0)=0, /(l) = ln|B|, B xy (l)=B xy , B xy (0) = 5{x - y), 



For calculation of \B\ we use the formula 

In \B\ = /(l) = f 1 /'(/.), /'(//) = SpS'O*) ■ (B.5) 
Jo 

Let us find the inverse matrix Note that the matrix B(fi) for < /i < 1 has the 

same analytic behaviour as the matrix B, therefore, for simplicity we only derive the inverse 
matrix B~ l . 

Show that matrix D 

D = L-E-J-E T -J (B.6) 
is inverse to B up to an exponentially small term 

B ■ D = 1 + 0(a r ) at r -> oo, 

where J XJ/ = 5(x + y — r), 

r dz x dz 2 (z 1 ) x (z 2 ) y r dz x dz 2 {z 1 )- x (z 2 ) v . . 

xy J\z lZ2 \<i 1 - z x z 2 P(zi 1 )Q(z 2 1 ) J\z 2 \<\ Zl \ Zi(z x - z 2 ) P(z 1 )Q(z 2 1 ) ' 

L* y = I = BL, if B(z^) = B~\z), 



., z B(z^ 

since 



and in our case B(z x ) = B 1 (z). 
At the first we calculate B ■ L 



(t> t\ I dzidz 2 B{zi) 1 - {z 1 /z 2 ) r+1 / 

^''wwW ri Zl ^Bw (zir _ t(z2)r _, (B8) 



|<1 2iZ 2j B(^2 ) 

Now one calculates the product B ■ E 

= I djzxz^r^z^ B( Zl ) l-( Zl /z 2 ) r+1 

^ '*v /|zi|<|2 3 |<|z2| ziz 2 (z 2 - z 3 ) P(z 2 )Q(z 3 l ) l-z 1 /z 2 

d(z 1 z 3 )(z 1 )-%z 3 )y B{z x ) , r d(z lZ3 )(zi)- x (z 3 )y Q(z^) 



zi\<\z 3 \ z x {z z -z x ) B{z 3 ) Jf|*i|<|*8| z l {z 1 -z 3 ) Q{z 3 x ) 

d(z 1 z 2 z 3 )(z 1 y-%z 3 )y(z 2 y +1 b( Zi ) 

a<\ Zi \<l (1 - Z\Z 2 ){\ - Z 2 Z 3 ) P(z 2 1 )Q{z 3 1 )' 



(B.9) 



In (|B.9|) the integration over z 2 yields a contribution which is propotional to a r+ and 
therefore, the last term in the right-hand side of (|B.9| ) is exponentially small for any x, y. 



The second term in the right-hand side of ( |B.9| ) gives the following contribution at the 
contraction of the contour over z 3 



d(z 1 z 3 )(z 1 )- x (z 3 )y Q(zT l ) 



Z3\<\zi\ 



Z\\Z\ - Z 3 



j— — 5(x — y) = -5(x-y), 



(B.10) 



where we took into account that the first term in the left-hand side of ( P-10|) reduces to zero 
at the extention of the contour over z 3 to one of infinite radius. 

Similarly for the first term in the right-hand side of ( P-9| ) one gets 



d(z 1 z 3 )(z 1 )- 3: (z 3 )y B( Zl ) 



+ 8{x - y). 



I\z 3 \<\zi\ z 1 (z 3 -z 1 ) B{z 3 ) 
Taking into account ( |B.10| ) and (|B.11|) , we obtain 

I d(z 1 z 2 )(z 1 )-%z 2 )y B( Zl ) 



(B.ll) 



>w J\z 2 \<\ Zl \ z 1 (z 2 -z 1 ) B(z 2 ) 

d{z 1 z 2 z 3 ){z 1 y-*{z 3 )y{z 2 y +l b( Zi ) 

\ Zi \<i {\ - z x z 2 ){\ - z 2 z 3 ) P(z 2 l )Q(z 3 l ) 

_ r d( Zl z 2 )( Zl y(z 2 )y B(z{ 1 ) 
zi-a/z! J\ zi z 2 \<\ (1 - z\z 2 ) B{z 2 ) 



f 



- 0(a r ) 



(B.12) 



From here it follows that B ■ E coincides up to 0(a r ) with the second term in the right-hand 
side of flBTgp of the product B ■ L. 



(B.13) 



Now let us consider the product B • JE T j 

B ■ JE T J = J ■ ( J ■ B ■ J ■ E T ) ■ J = J • (B T ■ E T ) ■ j, 

where we used J ■ B ■ J = B T . 

A transition from B, E to B T , E T means the following transformations in contour inte- 
grals (P-2| ) and (P-7|) : B{z) — > B(z~ l ), P{z) <-> Q{z) and therefore, no calculating, we have 
instead of ( gTg ) 



(B T ■ E 1 



xy 



d{ Zl z 2 ){ Zl Y{z 2 )y B{ Zl ) 



\ Zl z 2 \<i (1 - ziz 2 ) B(z 2 

d( Zl z 2 z 3 )( Zl y-*(z 3 )y(z 2 y +1 b(z^) 

/a<|«i|<l (1 - Z X Z 2 ){\ - Z 2 Z 3 ) Q(z 2 1 )P(z 3 1 ) 

and from (|B.13|) and (|B.14|) one gets 



B ■ JE T J 



d(z 1 Z2)(z 1 ) r - x (z 2 y- y B( Zl ) 

Bfzn 1 " 



0{a r ). 



[B.U) 



(B.15) 



'|ziz 2 |<i (1 — Zxz 2 ) ij y^ 2 

Collecting together Q, (1B~T2|) and (|B7T3|) , one proves that D = B' 1 + 0(a r ). Thus for 
f '((i) in (|B.5|) we obtain 

/'(//) = SpS'(/i) • - - J • E^/x) • J 



(B.16) 



Note that 

SpB' ■ J ■ E T ■ J = Sp J ■ B' ■ J ■ E T = Sp (B') T ■ E T = Sp {B') T ■ E T = SpB'(fx)T ■ £(//). 
From here one gets 



/V) = SpBV) (L(/i) - 2£(//)) . 

Let us denote by 

U x = Sp , C/ 2 = Sp (B'(fi)E(fi)) , /'(/i) = Z7i - 2C/ 2 . (B.17) 

At the first we calculate U\ 

dz\dz 2 



Ux = ^ <f —^—^z^ y+x B^(zi) \aB{z 1 )z 2 x+y B~> 1 (z 2 ) = 

x,y J Z \Z 2 

d Zl dz 2 B»( Zl ) \nB( Zl ) [1 - {z l /z 2 ) r+1 }{l - (z 2 / Zl Y +1 ] 
\zi\<\z 2 \ Z\Z 2 B»{z 2 ) {l-z l /z 2 ){l-z 2 /z l ) 

(r + l)^ln5(z)-(7 + 1 ^ ^ B ^ lnB ^ 



z \J\zi\<\z 2 \ J\^\<\zi\J (zi - z 2 ) 2 B^(z 2 ) 

dz x dz 2 B^ ( Zl ) In B ( Zl ) ( z^ 1 £ dz x dz 2 B^{z 1 )\nB{z x ) ( z 2 ^ r+] 



,\<\z 2 \ (z x - z 2 ) 2 B^(z 2 ) \z 2 J J122KI21I (zi - z 2 ) 2 B^(z 2 ) \z 1 
In this expression the last two terms are exponentially small (~ 0(a r ))) and therefore, 



Now let us calculate U 2 



f/ 2 = E/ 



d{z x z 2 z 3 _ y+x y x B^{z 1 )\nB{z x ) 



% I zi(l - z 2 z 3 y~ l ~ 2 3 P^z^Q^z, 1 ) ~ 

d( Zl z 2 z 3 B^ Zl )\nB{ Zl ) [1 - (zi* 2 ) r+1 ][l - (z 3 / Zl y +1 ] _ 
7|*8|<|*a|<M<i zi(l - z 2 z 3 ) P^(z 2 1 )Q^(z 3 1 ) (1 - ziz 2 ){\ - zs/zi) 

I d( Zl z 2 z 3 BxjzJlnBjz!) 

J\z 3 \<\z x \<\z 2 \ (zi ~ z 3 )(z 2 - zi)(z 2 ~ z 3 ) P^z^Q^ 1 ]) 
In the right-hand side of (|B.19| ) contracting the contour over z 2 , we obtain 

r d( Zl z 3 Q^z^)\nB{ Zl ) r d{z x z 3 B^ Zl )\nB( Zl ) 

j\*\<m (zi - z 3 y q^zs 1 }) Az S \<\ zl \ ( Zl - z 3 y b»(z 3 ) { h { ■ ) 

In the first term of the right-hand side of ( |B.20| ) let us move the contour over z 3 to 00. Then 



U 2 = jdzQ^z" 1 ) [in P{z) +\nQ(z~ 1 ) 



d_ { 1 

d~z \Q»{z-i) 



fi <f dz\nP(z)—\nQ(z- 1 ) + 0(a r 



(B.21) 



Substituting ( gig ) and ( ggg ) in flBT7| ), one gets 
f^) = (r + l)<f^lnB(z)-^ 



dz x dz 2 B^(z 1 )\nB(z 1 ) 



zi\<\z 2 \ J\Z2\<\ Z1 \) (zi - z 2 ) 2 B^(z 2 ) 

o 



2fi { dz\nP(z)—\nQ(z- 1 ) + 0(a r ). 
J dz 



(B.22) 



The difference of the integrals over z 2 in the second term of the right-hand side of ( |B.22| ) is 
a integral over z 2 around the point Zi, that is 



dz 2 







{z 1 ){z 1 -z 2 ) 2 B^{z 2 ) dz\B»{z) 



Z = Z\ 



Then the second term in the right-hand side of ( |B.22|) has the form 



H ldz\aB{z)—hxB(z) = - [ dz— fin 2 B(zj) = 0, 
J dz 2 J dz \ ' 

since the integral over a closed contour from a total derivative is equal to the zero. 
Thus 

f'(ji) = ( r + l) j — hxB{z)--2ii j> dz\nP(z)-^-\nQ(z- 1 ). 
Integrating this expression over // in the limits [0, 1], we obtain 

dz . _ , , r . . _ , , d 



(B.23) 



/(l) = ln|P| = (r + l) f-jhiB(z) - j> dz\nP{z) — lnQ^ 1 ). (B.24) 

Substituting here the explicit expressions for B(z), P(z) and from and (BjB) we 

have 

/ y In = v (j ^ In P(^) - / ^ In P(* -1 )) = i> (P(0) - P(0)) = 0, 

- Idz\nP(z)-^-\nQ(z- 1 ) = ^ 2 idz\nP(z) f — ?— ) = /r 2 (lnP(a) - lnP(/3)) . 

J dz " 7 V- 2 — a z — p J 



Thus 



(l-a 2 )(l-/3 2 ) 
(1 - a(3) 



(B.25) 
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